Determining Neutrino Mass from the CMB Alone 
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Distortions of Cosmic Microwave Background (CMB) temperature and polarization maps caused 
by gravitational lensing, observable with high angular resolution and high sensitivity, can be used 
to measure the neutrino mass. Assuming two massless species and one with mass m„ we forecast 
cr(rn v ) = 0.15 eV from the Planck satellite and a(rn v ) — 0.04 eV from observations with twice the 
angular resolution and ~ 20 times the sensitivity. A detection is likely at this higher sensitivity 
since the observation of atmospheric neutrino oscillations require AmJ > (0.04eV) 2 . 
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Introduction. Results from the WMAP Q] show the stan- 
dard cosmological model passing a highly stringent test. 
With this spectacular success of the CMB as a clean and 
powerful cosmological probe, and of the standard model 
as a phenomenological description of nature, it is timely 
to ask what can be done with yet higher resolution and 
higher sensitivity such as offered by the Planck instru- 
ments and beyond. In this Letter we mostly focus on neu- 
trino mass determination, with brief discussion of other 
applications. 

Eisenstein et al. |2( found that the Planck satellite can 
measure neutrino mass with an error of 0.26 eV. This sen- 
sitivity limit is related to the temperature at which the 
plasma recombines and the photons last scatter off of the 
free electrons, Td ec — 0.3 eV. Neutrinos with m„ < Td ec 
do not leave any imprint on the last-scattering surface 
that would distinguish them from m v = 0. 

Neutrinos with mass m v < Tdec would affect the ampli- 
tudes of gravitational potential peaks and valleys at in- 
termediate redshifts. Massive neutrinos can collapse into 
potential wells when they become non-relativistic, while 
massless ones freely stream out. The observed galaxy 
power spectrum (which is proportional to the potential 
power spectrum at sufficiently large scales), combined 
with CMB observations can be used to put constraints 
on m„ 3]. At present such an analysis yields an upper 
bound on m v of - 0.3 eV 0[H. 

The alteration of the gravitational potentials at late 
times changes the gravitational lensing of CMB photons 
as they traverse these potentials 0, Including the 
gravitational lensing effect, we find that the Planck error 
forecast improves to 0.15 eV. We also show that more 
ambitious CMB experiments can reduce this error to ~ 
0.04 eV. These mass ranges are interesting because the 
atmospheric neutrino oscillations require that at least one 
of the active neutrinos have m„ > 0.04 to 0.1 eV. More 
detailed considerations show that the sum of the active 
neutrino masses (which is what the CMB is most sensitive 
to) should be at least 0.06 eV. 

Tomographic observations of the galaxy shear due to 
gravitational lensing can achieve sensitivities to m v simi- 
lar to what we find here @, Q . Our work is distinguished 
by its sole reliance on CMB temperature and polarization 



maps which have different potential sources of systematic 
error. Complementary techniques are valuable since both 
of these will be very challenging measurements. 

The most stringent laboratory upper bound on neu- 
trino mass comes from tritium beta decay end-point ex- 
periments [l(ij | which limit the electron neutrino mass to 
< 2 eV. Proposed exeriments plan to reduce this limit 
by one to two orders of magnitude b y se arching for neu- 
trinoless double beta decay (/3P0v) [llj. A Dirac mass 
would elude this search, but theoretical prejudice favors 
(and the see-saw mechanism requires) Majorana masses. 
Like the CMB and galaxy shear observations, these fu- 
ture (3(3Qv experiments will be extremely challenging. 
Lensing of the CMB. The intensity and linear polariza- 
tion of the CMB are completely specified by the Stokes 
parameters, /, Q and U which are related to the unlensed 
Stokes parameters (denoted with a tilde) by X(n) = 
X(n + 8n) where X stands for /, Q or U. The deflec- 
tion angle, Sn, is the tangential gradient of the projected 
gravitational potential, 



n 



) = 2 / dr*(rn, r)(r - r s )/(rr s ) , 



(1) 



where r is the coordinate distance along our past light 
cone, s denotes the CMB last-scattering surface, n is 
the unit vector in the n direction and * is the three- 
dimensional gravitational potential. 

The statistical properties of the /, Q and U maps are 
most simply described in the transform space: cltQ), 
<Ze(1), and ds(l) where ot is the spherical harmonic 
transform of / and cie and olb are the curl-free and 
gradient-free decompositions, respectively, of Q and U 
Jl2l Il3| . In this transform space the effect of lensing 
by mode 4>(L) (harmonic transform of 0(n)) is to shift 
power from, e.g., ot(L — 1) to aril)- Lensing also mixes 
cle into as and any as into cle |lj], thus generating 
scalar B (curl) mode correlations. 

Lensing smoothes out the features in the two-point 
functions, also called angular power spectra, C" a , where 
{a a (l)a* a ,{Y)) =C? a ' 5{\-\')/[2ttI{1 + \)] and a stands for 
T, E, or B |5j. As explained later, in our analysis we use 
the unlensed power spectra, C" Q . The information from 
lensing is added through the two-point function of the 
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lensing potential, (0(L)0*(L')) = Cf^5(L-L , )/[27^L(L + 
l)], which can be inferred from the temperature and po- 
larization map 4-point functions |15| . In Figurc^we plot 
the deflection angle power spectrum, Cf d = 1(1 + 1)C, . 
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FIG. 1: Top panel: Deflection angle power spectrum Cf d for 
the fiducial model (m v = 0). Bottom panel: 100xdCf d / dm v x 
(Am„/Cf d ) (dark) and 100 x dCf d /dw x x (Aw x /Cf d ) (light) 
for Am„ = 0.1 eV and Aw x = 0.2. 



We calculate the 2-point functions using a publicly 
available code, CMBfast 0, which was modified to in- 
clude a scalar field dark energy component, to calculate 
Cf d , and to include the effect of massive neutrinos on the 
recombination history (through the expansion rate). We 
use the Peacock and Dodds prescri ptio n to calculate the 
non-linear matter power spectrum 

Effect of neutrinos. The lower panel in Figure ^ shows 
the differences in the power spectra between our fiducial 
model and the exact same model but with one of the 
three neutrino masses altered from zero to 0.1 eV. The 
error boxes are those for CMBpol (described below; see 
Table 1). The Cf d are noise-dominated at I > 600 for 
CMBpol. 

The signature of a 0.1 eV neutrino in the angular power 
spectra, in the absence of lensing, is at the 0.1% level. 
Such small masses are only detectable through their ef- 
fect on lensing, which comes through their influence on 
the gravitational potential. Replacing a massless compo- 
nent with a massive one increases the energy density and 
therefore the expansion rate, suppressing growth. The 
net suppression of the power spectrum is scale depen- 
dent and the relevant length scale is the Jeans length for 
neutrinos 0, 0, which decreases with time as the 



neutrino thermal velocity decreases. This suppression 
of growth is ameliorated at scales larger than the Jeans 
length at matter-radiation equality, where the neutri- 
nos can cluster. Neutrinos never cluster at scales smaller 
than the Jeans length today. The net result is no effect 
on large scales and a suppression of power on small scales, 
resulting in the shape of 5Cf d /Cf d in Figure ^ 
Error forecasting method. 

The power spectra we include in our analysis are Cf T , 
Cf E , Cf E (unlensed), and Cf d . We do not use the lensed 
power spectra to avoid the complication of the correlation 
in their errors between different £ values and with the er- 
ror in Cf d . Using the lensed spectra and neglecting these 
correlations can lead to overly optimistic forecasts [2fl |. 
If we include the lensed spectra instead of the unlensed 
ones, the expected errors on w x and m„ for CMBpol (see 
Table 1) shrink by about 40% and 30% respectively. 

The distortions to the angular power spectra due to 
a 0.1 eV neutrino and changes of order 10% in w x are 
very small. We have taken care to accurately forecast 
the constraints possible in this mass range. First, we 
make a Taylor expansion of the power spectra to first 
order in all the cosmological parameters. Then, given the 
the expected experimental errors on the power spectra, 
the expected parameter error covariance matrix is easily 
calculated. 

The Taylor expansion works better and suscepti- 
bility to numerical error is reduced with a careful 
choice of the parameters used to span a given model 
space 0, EH S E3- We take our set to be V = 



{Uj m , UJb,UJ,y, 



, fc 3 P|(fc/),n s ,n' s , j/hc}, with the 



assumption a flat universe. The first three of these are 
the densities today (in units of 1.88 x 10 _29 g/cm 3 ) of cold 
dark matter plus baryons, baryons and massive neutri- 
nos. Next two are the angular size subtended by the 
sound horizon on the last-scattering surface and the ra- 
tio of dark energy pressure to density. The Thompson 
scattering optical depth for CMB photons, r, is parame- 
terized by the redshift of reionzation z^. The primordial 
potential power spectrum is assumed to be fc 3 P|(fc) = 
k 3 f Pi(k f )(k/k s ) n °- 1+n 's ln{ - k / k ^ with kf = 0.05 Mpc" 1 . 
The fraction of baryonic mass in Helium is j/hg • We Tay- 
lor expand about V = {0.146, 0.021, 0, 0.6, -1, 6.3, 6.4 x 
10" 11 , 1,0, 0.24}. 

We follow [24| to calculate the errors expected in Cf T , 
CT E and C EE given Table 1 . For errors on Cf d we follow 
hj. The errors on the unlensed spectra in the regime 
where lensing is important (deep in the damping tail) 
are certainly underestimated because reconstruction of 
the unlensed map from the lensed map will add to the 
errors. However, this is not worrisome since we limit all 
the unlensed spectra to I < 2000, and a further restric- 
tion to I < 1500 (where lensing is least important) only 
increases the error on m„ by about 10% for CMBpol. 
Experiments. We consider Planck |25j . a high-resolution 
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Experiment Sax fi£ v (GHz) 9 b A T A P 



Planck 


2000 


2500 


100 


9.2' 


5.5 


CO 








143 


7.1 
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11 








217 


5.0' 


13 


27 


SPTpol (/ sky = 


: 0.1) 2000 


2500 


217 


0.9' 


12 


17 


CMBpol 


2000 


2500 


217 


3.0' 


1 


1.4 



TABLE I: Experimental specifications. We use the unlensed 
spectra (Cf T , Cf E , Cf E ) only at I < 2000. For (j> reconstruc- 
tion we use only data with I < lmax.' B ' ■ 

version of CMBpol 01 > an d a polarized bolometer array- 
on the South Pole Telescope |42| we will call SPTpol. 
Their specifications are given in Table 1. We assume 
that other frequency channels of Planck and CMBpol 
(not shown in the table) will clean out non-CMB sources 
of radiation perfectly. Detailed studies have shown fore- 
ground degradation of the results expected from Planck 
to be mild [H [H HI . At 1 > 3000 emission from dusty 
galaxies will be a significant source of contamination. 
The effect is expected to be more severe for temperature 
maps. Hence we restrict temperature data to I < 2000 
and polarization data to I < 2500. 

Results. We emphasize the ability of the experiments to 
simultaneously determine PJ,, w x an d m v E3- These 
all affect the amplitude of P$ at late times, the latter 
two due to their effect on the rate of growth of density 
perturbations. If we were only sensitive to the amplitude 
of Cf d then there would be an exact degeneracy between 
these three parameters. However, the Z-dependence of 
the response of Cf d to these parameter variations breaks 
this would-be degeneracy, allowing for their simultaneous 
determination. 

The effect of m„ can easily be disentangled from that 
of w x . We have already discussed the /-dependence of 
d\nCf d jdm v shown in Fig. ^ as resulting from the 
scale- and time-dependence of d In P$,/dm u . The l- 
dependence of d\nCf d / dw x has the opposite sense. Al- 
though the suppression of P$ for increasing w x is nearly 
/c-independent, the effect is larger at late times — hence 
the radial projection gives a larger effect at low I. The 
effects of m„ and w x are sufficiently distinct to allow for 
their simultaneous determination. We point out that the 
effect of w x is more pronounced for larger values due to 
two reasons. One, dark energy starts to dominate earlier 
(which implies larger uniform suppression) and two, per- 
turbations in dark energy on large scales are enhanced 
for large w x . 

The difference in the response of Cf d to m„ and w x 
allows for, e.g., Planck to detect the acceleration of the 
Universe (w x < — l/(30 x )) at the 2a level. Such a con- 
firmation would be valuable given the deep theoretical 



implications of acceleration 30]. Hu 20] has previously 
noted this result obtained with the assumption m„ = 0. 

As is well known, the P| can be determined indepen- 
dently of the lensing signal, through use of a signal at 
large angular scales. One combines C EE and Cj E at 
I < 20 where t hey are proportional to PJ,t 2 and PJ>t 
respectively jHHa] with the TT, EE and ET spectra at 
20 < I < 2000 where they are proportional to P^e~ 2r . 

If we assume a single-step transition for the ioniza- 
tion history Planck can achieve ct(t ) = 0.005 2]. How- 
ever, foreground contamination |27|. and modeling un- 
certainty in the ionization history [3l| can increase this 
uncertainty. For these reasons we conservatively ignore 
polarization data at I < 30 and instead set a prior, by 
hand, of <t(t) = 0.009; including the I < 30 polarization 
data would (perhaps artificially) achieve a smaller <j(t). 
In the end, r is determined (only slightly) better than 
this prior because there is some constraint on PJ> from 
the lensing signal. Note that since Pf,e~ 2T is so well- 
determined, we always expect <r(lnPj>) = 2<j(t), as we 
find.. 

An extended period of reionization, as suggested by 
the combination of WM AP and quasar observations [32j , 
may have large spatial fluctuations in the ionization frac- 
tion. Such "patchy" reionization would lead to a large 
diffuse kinetic SZ contribution to Cf T at high I 013, 
possibly larger than the lensing contribution. Fortu- 
nately the analogous effect in the polarization is much 
smaller. For a conservative upper bound on how patchy 
reionization could degrade er(m„), we restrict the tem- 
perature data to I < 1000 and find a(m v ) = 0.045 eV for 
CMBpol and 0.34 eV for Planck. 

The primary motivation for CMBpol is the detection 
of the B mode due to gravity waves produced in infla- 
tion. The amplitude of this signal would directly tell 
us the energy density during inflation. Following the 
calculation in [H |H we find a 3<7 detection is possi- 
ble for CMBpol if the energy density during inflation is 
greater than p min = (2 x 10 15 GeV) 4 ; p^ n is an order of 
magnitude smaller than the GUT scale. We note that 
Pmin oc 1/r, approximately, for 0.05 < r < 0.2 and we 
have assumed r = 0.1. This scaling with r suggests that 
the reionization feature in the B mode at the largest an- 
gular scales is important and therefore a full-sky experi- 
ment is necessary to achieve this sensitivity level. 

The scalar spectrum determined from high-resolution 
CMB observations (the constraining power comes from 
primary CMB) can also be a useful probe of inflation, as 
studied recently by [33 ■ If ns — 1 = 0.07, the central 
value in fits to WMAP and other observations 0, then 
inflationary models generically predict n' s ~ (ns — l) 2 = 
0.005 which will be detectable at the 3c level by CMBpol. 

Determining u>b and yue to high precision will facilitate 
precision consistency tests with Big Bang Nucleosynthe- 
sis (BBN) predictions. It will also be useful in constrain- 
ing non-standard BBN. For example, determining u>b and 
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TABLE II: 
Error Forecasts 







Experiment 


ra v (eV) 


w x 


In P|, 


ns 


n' s 


9s (deg) 


T 


In ui m 


In uib 


Z/Ho 


Planck 
SPTpol 
CMBpol 


0.15 
0.18 
0.044 


0.31 
0.49 
0.18 


0.017 
0.018 
0.017 


0.0071 
0.01 
0.0029 


0.0032 
0.006 
0.0017 


0.002 
0.0026 
0.00064 


0.0088 
0.0088 
0.0085 


0.0066 
0.0087 
0.0022 


0.0075 

0.01 
0.0028 


0.012 
0.017 
0.0048 







NOTES.— Standard deviations expected from Planck, SPTpol and CMBpol. 



yue to high precision allows strong constraints to be put 
on the number of relativistic species N (or equivalently 
the expansion rate) during BBN. If a(yH e ) is small, then 
a(N) = cr(y ffe )/0.013, which for CMBpol works out to 
cr(iV) = 0.4. Constraints on N have important repercus- 
sions for neutrino mixing in the early universe, and hence 
on neutrino mass models [ssj . 

Conclusions. Gravitational lensing of the CMB is a 
promising probe of the growth of structure and the fun- 
damental physics that affects it. High sensitivity, high 
resolution maps will allow us to measure the lensing sig- 
nature well enough to simultanteously constrain m„, w x 
and P$. A future all-sky polarized CMB mission aimed 
at detecting gravitational waves is likely to succeed in 
determining neutrino mass as well. 

We thank J. Bock, S. Church, W. Hu, M. 
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useful conversations. 
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